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We present a simple, analytic model for pump-probe spectroscopy in dilute atomic gases. Our 
model treats multilevel atoms, takes several broadening mechanisms into account and, with no free 
parameters, shows excellent agreement with experimentally observed spectra. 

PACS numbers: 32.30.-r 32.70.Fw 



I. INTRODUCTION 



There has been much interest recently in the effect of 
the propagation of a laser beam through dilute vapors of 
alkali atoms [J, [3, H, 0, 0] due to their application in laser 
cooling |1, 0| , chip-scale atomic clocks and magnetome- 
ters P, 'QI, laser stabilization [l^, fill. fT3 |. and electro- 
magnetically induced transparency In order to un- 
derstand the spectra, several models have been proposed 
which accurately predict the Doppler-broadened and sub- 
Doppler absorption spectra of a dilute vapor with weak 
and strong probe beams [H, Hi, [S [3 [13, 111 • The com- 
mon approach of numerically solving the semi-classical 
density matrix equations, however, requires intensive 
computation This may be unnecessarily compli- 

cated for several applications in which coherence effects 
are negligible and beam powers relatively l ow, such as 
studies of optical pumping ^ j transit- [l|, [3, [i3| and 
power-broadening and number density measurements. In 
this paper we show that a simple ab-initio model based 
upon rate equations can predict both Doppler-broadened 
and sub-Doppler spectra with high accuracy for a large 
range of beam powers and widths. 

Our model has been developed primarily to focus upon 
dilute vapors of alkali atoms, specifically rubidium, but 
the results should be generally applicable to any dilute 
vapor. We use the Einstein rate equations to calculate 
the steady state population densities and therefore spec- 
tra are assumed to be measured on timescales greater 
than any coherence effects [Tsl]. The model also as- 
sumes coUimated, pump and probe beams, whose spec- 
tral widths are less than the spontaneous decay rates of 
the atoms. 

This paper begins with a short review of the theory 
behind saturated absorption spectroscopy; the next sec- 
tion covers the effects of optical pumping on absorption 
lineshapes with multilevel atoms, which accounts for the 
majority of observed features. We then compare the re- 
sult of this simple model with experimental pump-probe 
spectra of rubidium and show that the fit is accurate, 
even without any free parameters. 



II. PUMP-PROBE THEORY 

The absorption of a weak probe beam propagating in 
the z direction through a dilute gas is characterized by 
the Beer-Lambert relation Il9l. 
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where Ny is the number density of atoms, Iq is the in- 
cident probe intensity and ct(w) is the absorption cross 
section for the electric dipole transition z ^ fc, in which 
i is the ground state and k is the excited state, and the 
difference in state energies E^. — Ei = Hluh^: 
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where Ui and Uk are the fractional populations of the 
ground and excited states, Bik and Bki are the Einstein 
coefficients for the absorption and stimulated emission of 
a photon, respectively, and 
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is the normalized Lorentzian function which character- 
izes the atom's response to the incident field, where lo 
is the frequency of the laser and 7 = 1/r is the natural 
linewidth of the transition which has the excited state 
lifetime r. 

In common with many texts, we set the excitation rate 
to be proportional to the energy density in each spectral 
mode p{uj) for a broadband light source. For high reso- 
lution spectroscopy we require the measuring tool — the 
laser — to have a finer resolution than the subject under 
investigation and therefore the spectral linewidth of the 
laser must be Acj < l/r. Hence the spectral energy den- 
sity of a nearly monochromatic beam with an electric 
field amplitude S may be written as [l^] 



p{u)^\e£^L{u) = -L{Lo), 
I c 



(4) 



which has already been assumed in Equation [51 For 
states with degeneracy gi^k, the Einstein coefficients are 
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where e is the electron charge, eo is the vacuum permit- 
tivity and luikl = \fiki\ = Cik\{Jk\\r\\Ji)\ are the dipole 
matrix elements for electric dipole transitions between 
states with spin-orbit angular momentum Ji^k, in which 



|(^fe||r|| Jj)| = flo 
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is the reduced dipole moment [1, [2l| and Cik are the 
Clebsch-Gordan coefficients for the mp sublevels which 
define the relative strength for each transition. For unpo- 
larized beams with no external magnetic fields (no quan- 
tization axis) the dipole matrix elements are; 



(7) 



Values of J2mp ^fk '^^^ ^e obtained from tabulated re- 
sults 22] or via a lengthy calculation [23] . par It is com- 
mon in many texts to include a prefactor of 1/3 in Equa- 
tion O to average over the 3 possible orientations of the 
atom with no quantization axis [l^, [13, |25]. For degen- 
erate states, the summation over all mp sublevels yields 
the same value for all polarizations and therefore in Equa- 
tion 3] we assume linear polarization and avoid the 1/3 
prefactor. 

The lineshapes of the resonances will be inhomoge- 
neously broadened by the motion of the atoms. For a di- 
lute vapor, the atoms have a mean free path greater than 
the dimensions of the cell, therefore collisions are negligi- 
ble and the spread of velocities will follow a Maxwell- 
Boltzmann distribution with a mean temperature, T, 
equal to that of the cell walls. We therefore have the 
distribution of atoms with velocity v, 
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and mean speed, 
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where M is the atomic mass and fcs is the Boltzmann 
constant. The change in intensity for a weak probe beam 
through a dilute vapor is therefore 

d/(^) huj,k7 ^. , „ /+°" fDiv){n^~fnk) 

- NyloBik I rdw, 



dz 
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(10) 

where kv = {uj/c)v is the Doppler shift which has been 
included into the Lorentzian lineshape 
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and similarly with p(lu,v). The integrand has a Voigt 
lineshape which cannot be solved analytically but may 
be calculated from tabulated values Q . 



A. Saturation 

The common experimental arrangement is for a strong 
pump beam and weak probe beam, derived from the 
same laser and hence equal in frequency, to propagate 
in opposite directions through the atomic sample. The 
pump beam affects the population difference, rii — ^nk, 
and thus reduces the probe absorption near resonance. 
Only atoms which travel nearly perpendicular to both 
beams, and thus have a zero velocity component along 
the beam direction, will be pumped and probed simul- 
taneously, and therefore narrow 'Lamb dips' form in the 
Doppler-broadened profile at resonance |25| . 

Using the Einstein rate equations for a closed, degen- 
erate two level system (n; + nk = I, gk ~ 9%), the steady 
state population difference is 



m - Uk 



2BkM^, -v)If 
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where Ip is the pump beam intensity, and 



(12) 
(13) 

(14) 



is the saturation intensity at which nk = 1/4 on res- 
onance and the negative sign for the velocity in Equa- 
tion [13 reflects the counter-propagating geometry. As 
the pump intensity increases, the population difference 
tends to zero, the atomic sample becomes transparent to 
further excitation, and thus the probe beam shows de- 
creased absorption. 

This model is commonly presented in text books 
(T^ [2^ and is known as saturated absorption spec- 
troscopy. Unfortunately it fails to model realistic systems 
in two important ways. Firstly, the system is not closed 
as atoms may enter or leave the beams; this limits the 
interaction time and re-thermalizes the population dif- 
ference. Secondly, closed two-level systems are very rare 
and the atoms will generally have several ground and ex- 
cited states into which they may be optically pumped. 
If these additional states are separated by frequencies 
greater than the laser linewidth, then once these states 
are populated the atom will become transparent to the 
laser radiation. This process can occur even at low pump 
intensities and therefore is a significant factor affecting 
the depth of the sub-Doppler features [1] . 



B. Optical Pumping 

The simplest multilevel system involves a single excited 
state ]c) which may decay into two ground states \a) and 
\b), as shown in Figured] The excited state decays at 
a total rate 7c but the fractional rates into each ground 
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FIG. 1: (colour online) Simple 3 level system with optical 
pumping on \a) \c) transition. The steady state popula- 
tions in each state depend upon the decay fractions and 
the populations of atoms entering and leaving the beam at rate 

ro. 



state (branching ratio) depend upon the ratio of Clebsch- 
Gordan coefRcients, 7ca and 'ycb, where 



C, 



^ac + ^bc 



(^ac + C^bc 



7c- 
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Atoms drift into and out of the beam at a rate Fq, 
entering with the equihbrium population distribution and 
leaving with the optically-pumped distribution. This is 
discussed further in Section Hi Dl 

We consider the absorption from state |a), so that any 
population excited from \b) is negligible. The rates of 
change of the state populations are therefore 



dt 



+nc{Bcapi0J, V) + 7ca) + Afa^O (17) 

^ = ro(-nfa+AAfc)+n,7,fc (18) 
dt 

duc 

-nciBcapioJ,v)+-fc + To), (19) 
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where A/'a and A/{, are the fractional equilibrium popula- 
tions of states \a) and respectively, when Bac — 0. 

The equilibrium populations depend upon the thermal 
distribution of energies [10, |13l , but for most atoms the 
ground state splitting is of the order GHz, while the ex- 
cited state is separated by hundreds of THz from the 
ground state. The mean energy per degree of freedom is 
^ksT, which at room temperature corresponds to a fre- 
quency of THz and so we may assume the ground states 
to be equally populated and the excited state population 
to be negligible. The equilibrium population of state \a) 
is therefore 
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When the rate of the laser scanning across the resonance 
is much slower than any depopulation mechanisms in 
the system, we may assume that the populations have 
reached their steady state values (ria — hh = he = 0) and 
hence become 

nc{BcaP{'^,v) +lca) +J^a^O /^l ^ 

rin = ; ^ , (21) 
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where we assume |w — ujca\ <^ \uj — LUcb\ and hence may 
neglect n^. This makes the derivation valid for any open 
two-level system. We then substitute Equation [22] into 
inland rearrange to find ria, 
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is the saturation parameter, and 
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is the optical pumping parameter which enhances the 
saturation of the spectra. As derived in Section III Al the 
absorption cross section for the incident beam is propor- 
tional to the population difference between the ground 
and excited states, BacTia — Bca^c which we may find in 
terms of Afa by using equation [22] 



Bacna — Bcaflc — BacMa 



1 
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BacAfa^N 
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where AA^ is the population difference caused by the 
pump beams. This may be rearranged into a form similar 
to Equation [T51 
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in which is the reduced saturation intensity reported 
in the literature (T^, [l^l which is valid for an open two 
level system: 



— T 



1 + Tohk 



2 + (7/c - 7/ci)/ro 
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The numerator in the first term of Equation 1271 accounts 
for the steady state population pumped into the excited 
states and the denominator accounts for the population 
pumped into the ground states. 
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the above model are given by 



FIG. 2: (colour online) The effect of optical pumping in 
multilevel atoms. The dashed line shows the population dif- 
ference using Equation \27\ i.e., no optical pumping. The 
solid blue line shows the population difference for the **^7?6 
_F = 3 ^ -F' = 3 transition in which optical pumping sig- 
nificantly empties the upper ground state at intensities much 
lower than saturation. The solid red line shows the popula- 
tion difference for the *^7?6 F = 3 — > F' = 4 transition which 
is normally considered 'closed'; we would expect the popula- 
tion to follow the dashed line and equal 1/2 at the satura- 
tion intensity. In the absence of off-resonant pumping the 
population will tend toward an equilibrium with the transit 
rate To (AN ~ 0.44 in this instance) and so the population 
at saturation would equal one half of this value from unity 
(1 — |0.44 = 0.72^. The equilibrium value is not reached at 
high intensities due to increased off resonant optical pumping, 
but evidence for this can be seen in the slight shoulder around 
Ip/Uk = 100. 



C. Multilevel atoms 



Equation [26] may be intuitively extended to include 
multiple ground (i) and excited (fc) levels. 



7' [ Ip ll"^ 
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and 



7/2 



(33) 



(34) 



Note that for Equation[33]the Doppler shift for the probe 
beam is opposite to that of the pump beam. 

This model may be simplified by considering the steady 
state population for open or closed transitions with fast 
(MHz) decay rates. Off-resonant pumping into the dark 
state must be included both to model crossover peaks cor- 
rectly and because an equilibrium is set up between 
pumping and transit of atoms into the beam (see Figure 
[2]). For an open transition the atoms will be pumped into 
the dark state before a significant population builds up 
in the excited state, so one may use the following form 



A7V 



(35) 



For a closed state (such as ^^Rb F = 3 — > F' = 4) a sig- 
nificant population will build up in the excited state, but 
the off-resonant excited state populations will be negli- 
gible. In such a situation off-resonant optical pumping 
into the dark ground states reduces the equilibrium pop- 
ulation as shown in Figure [2] and Equation 1291 mav then 
be simplified to the form 



AiV : 



1 - OLi, 



(36) 
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(30) 
(31) 



The absorption of the probe beam is therefore 



D. Transit broadening 

Transit broadening is normally considered in the con- 
text of a limit to the interaction time and hence to the 
fundamental resolution of the spectra. For a Gaussian 
laser beam with a l/e^ radius R passing through a dilute 
gas of atoms with a mean velocity u, the standard result 
[23 , [26j for the transit rate of an atom through the beam 
is 



^ = -^Ny fyo{v)ANi-v)J2B^,piu;,v)dv. 

(32) 

Because of optical pumping even at low intensities, the 
pump beam can no longer be assumed to have a negligible 
effect upon the spectra. Therefore, the pumping rates in 
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However, we can see in Equation]^ that the transit rate 
has a significant effect on the optical pumping rate. It is 
common practice to use large beams to reduce the tran- 
sit broadening but this has the effect of increasing the 
time during which atoms may be optically pumped into 
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the dark state, and hence acts to decrease the resohition. 
In Figure [2] we see that optical pumping occurs at in- 
tensities much lower than 7^^; in fact, the atom can be 
pumped into the dark state \b) with the absorption of a 
single photon and so this complicates the calculation of 
the beam width. If we assume that significant optical 
pumping occurs at a specific intensity, and calculate 
the width of beam at this intensity as the total pump 
intensity is changed, the transit rate becomes; 



the counter-propagating beams cross at a non-zero angle 
6. This results in a broadening [2^ 



(fcp -ko) -v 



2uuji 



sin(i( 



(39) 



where kp and ka are the pump and probe wavevectors, 
respectively. 



- To 



if Ip > 
if Ip < 



(38) 



The transit rate here depends upon the specific transition 
{i —^ k), since the apparent increase of beam width will 
depend upon the individual transition strength. It is not 
obvious what value of I^p of should be used; however in 
the following model we find empirically that the reduced 



saturation intensity (/^ 



fits the data well. 



E. Additional broadening mechanisms 

The experimental setup may introduce extra broaden- 
ing of the spectral features and although they are usually 
much less than the natural and Doppler linewidths, they 
can be on the order of the transit rate. 

The following broadening mechanisms have a 
Lorentzian lineshape or effect the frequency detuning, 
therefore may be included to our model by summation 
with the natural decay rate [20| . 



3. Beam collimation 

An uncollimated beam passing through the vapor cell 
will result in a variation of transit time along the cell. 
Wavefront curvature also broadens the Lamb dip in much 
the same manner as geometrical broadening. To limit 
this effect the radius of curvature of the wavefronts must 
be much greater than kR^/A [25|. 



4- Collisional broadening 

We have limited our model to dilute gases in which 
collisions during the interaction time are negligible. In 
dense gases the collisional cross-section can be as large 
as the absorption cross-section and can affect the spectra 
by broadening and shifting the peak. Collisions may also 
affect the distribution of population amongst the ground 
states. The magnitude and nature of the effect depend 
upon density, temperature and collisional partners (i.e., 
non-identical atoms in the case of a buffer gas) [l^l • 



Laser linewidth 



III. RUBIDIUM 



As mentioned earlier, the resolution of the spectral fea- 
tures is constrained by the tool wc use to measure them: 
the laser. Diode lasers are a common tool for pump probe 
spectroscopy and may have linewidths down to the hun- 
dreds of kHz region with the addition of an external cav- 
ity. We assume a Lorentzian spectral laser line-shape 
with a full width at half maximum (FWHM) Fl in an- 
gular frequency units. 



2. Geometrical broadening 

The sub-Doppler resolution depends upon atoms 
traversing the beams at a perpendicular angle. There- 
fore, if the beams are not exactly anti-parallel they will 
sample a non-zero atomic velocity component and hence 
decrease the resolution. In the experimental arrangement 
described in Section HVl perfect beam overlap is achieved 
using a polarizing beam-splitting cube. For some experi- 
ments the polarizations of each beam may need to be con- 
trolled independently, and back-reflections into the laser 
cavity may be unwanted, favoring a geometry in which 



Natural rubidium occurs in two stable isotopes, ^^Rb 
and ®^Rb, with fractional abundances, //, of 0.2783 and 
0.7217, and atomic masses 84.912 and 86.909, respec- 
tively 27]. The ground state is 5Si/2, with nuclear mo- 
ment 7 = 5/2 for s^Rb and 7 = 3/2 for s'^Rb, and 
orbital angular momentum L — and spin S = 1/2. 
The excited state under investigation is 5P3/2 state with 
7 = 1 and the transition between these states is known 
as the 7?2 line. Due to the non-zero nuclear angular mo- 
mentum each state is split into 2 J -I- 1 hyperfine levels, 
each of which has 2(27 -1-1) magnetic sublevels. Fig- 
ure [3] shows the hyperfine structure of both isotopes and 
the energy splitting between levels. The Clebsch-Gordan 
coefficients, saturation intensity (Equation I14p and frac- 
tional decay rates "fik used in the model are shown in 
Table n The natural lifetime is 26.25 ± 0.07 ns [H] for 
both isotopes. 

Rubidium has a meltingpoint of 39.3°Cand the number 
density is then given by [29| 



ksT 



133.323 X IQP^ 



(40) 
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FIG. 3: Hyperfine structure of^^Rb and ^'^ Rh. 
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FIG. 4: The pump-probe spectroscopy apparatus used for the 
spectra in SectionW[ Approximately 3mW is picked off from 
the external cavity diode laser (EC'DL) with a beam-splitter 
plate (BSP). The ECDL is protected from back reflection by a 
Faraday Optical Isolator (FOI), so that the beam may be per- 
fectly overlapped using a polarizing beam splitter cube (PBSC). 
The linear polarization of the probe (left to right through the 
vapor cell) and pump (right to left) beams are rotated by half 
waveplates (HWP) so that the pump is transmitted and the 
probe reflected by the PBSC. The power of the pump beam is 
controlled with a variable neutral density filter (VND). 
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TABLE I: Rubidium model parameters for linear polarization 
and degenerate hyperfine levels. 



where the exponent subscript j — S, L corresponds to 
sohd and liquid rubidium, respectively, with 



PS 



-94.04826- 



1961.258 
f 



-0.03771678 x T + 42.57526 ln(r) (41) 
4529.535 

PL = 15.88253 

+0.00058663 x T - 2.99138 ln(T) . (42) 



in the Littrow configuration with a beam diameter of 
2.1 ± 0.1 mm in the vertical plane. The beam was ellipti- 
cal with a horizontal width equal to twice the height; 
the greatest transit broadening thus results from the 
narrower dimension and the model therefore uses the 
above value. The beams were overlapped through the 
vapor cell by sending the pump beam through a polariz- 
ing beam splitter cube (PBSC); the counterpropagating 
probe beam had a linear polarization perpendicular to 
the pump beam and was reflected by the PBSC onto the 
photodiode. This layout allowed perfect overlap of the 
beams and a reduced footprint of the apparatus. 

The laser frequency was scanned via rotation of the ex- 
ternal cavity grating with a piezoelectric transducer and 
spectra were averaged over three scans. The scans were 
not linear, so the frequency axis was calibrated using the 
12 resonance and crossover peaks (see Figure [3]) by fit- 
ting to a fourth-order polynomial; the standard deviation 
of the fitted peaks from tabulated values was < 400 kHz. 
The probe power was kept at 1/iW and the laser linewidth 
was 1.0 ± 0.2 MHz [28]. The vapor cell had length of 
75 mm and a temperature of 18.0 ± 0.1 °C. No attempt 
was made to null the Earth's or nearby magnetic fields. 



COMPARISON BETWEEN EXPERIMENT 
AND THEORY 



IV. PUMP-PROBE APPARATUS 

We have tested our model experimentally with rubid- 
ium vapor in a pump- probe apparatus, normally used to 
stabilize diode lasers for a magneto-optical trap, shown 
in Figure 4. The diode laser had an external cavity 



The experimentally measured pump probe spectra for 
the upper ground states of ^^Rb and *^Rb are shown to- 
gether with predictions of the theoretical model in Figure 
|6li, where i=a,b,c,d relate to the pump powers 20 /iW, 
100 AiW, 500 /.tW and 2500 /xW, respectively. Figure [H 
is split into four parts: plots i-3 shows the full absorp- 
tion spectrum with the experimental data in black and 
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FIG. 5; (colour online) The effect of beam width on sub- 
Doppler features. The dashed line and solid red line corre- 
spond to the Figure\^-2. The blue line shows the same m,odel 
but without the additional beam broadening factor from Equa- 
tion \S8\ The green line includes the additional broadening 
factor but assumes a beam width twice that in the red 

plot. The laser power remains the same for each plot. 

the theoretical curves in red, plots i-1,2 show a magni- 
fied sections of the Lamb dips (i-l=^^Rb, i-2=^^Rb) and 
plots z-4 shows the experimental data subtracted from 
the theoretical curves (residuals). 

It is apparent that our theoretical model accurately 
predicts the height and width of each absorption line. 
The residuals have a standard deviation of less than 1%, 
and much of this is due to experimental noise and cal- 
ibration errors. In Figures 6a-4 and 6b-4, the noise is 
dominated by a sinusoidal signal from power line pickup. 

Our model appears to slightly overestimate the height 
of the Lamb dips at low powers. This may be due to 
errors in the beam power measurement or the elliptical 



cross section of the laser beams. This latter point compli- 
cates the effect of beam width on optical pumping, and 
the arbitrary value used for I^p in Section lll Dl The effect 
of these values is shown in Figure [5l 

VI. CONCLUSION 

A simple model of pump-probe spectroscopy based 
on the rate equations has been presented and compared 
against the experimental spectrum of rubidium. Our 
model can describe multilevel atoms and fits the data 
well with the residual difference less than 1% over a large 
range of pump powers. The model is valid for any dilute 
gas when coherent effects are negligible, and accounts for 
finite laser linewidth, optical pumping and transit time 
broadening. 

The most significant effect upon the spectral features 
is the optical pumping during the atom's transit across 
the beam. As opposed to saturation broadening, which 
is due to the equalization of populations in the ground 
and excited states preventing further absorption, opti- 
cal pumping can transfer population into dark states, 
which also prevents further absorption but may occur for 
a single photon absorption. The time taken by the atom 
to transverse the beam significantly affects the optical 
pumping and therefore careful attention must be made 
in defining the beam width. 
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